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We present cross section measurements for Z/7*+jets+X production, differential in the trans- 
verse momenta of the three leading jets. The data sample was collected with the DO detector at the 
Fermilab Tevatron pp collider at a center-of-mass energy of 1.96 TeV and corresponds to an inte- 
grated luminosity of 1 fb~^. Leading and next-to-leading order perturbative QCD predictions are 
compared with the measurements, and agreement is found within the theoretical and experimental 
uncertainties. We also make comparisons with the predictions of four event generators. Two parton- 
shower-based generators show significant shape and normalization differences with respect to the 
data. In contrast, two generators combining tree-level matrix elements with a parton shower give a 
reasonable description of the the shapes observed in data, but the predicted normalizations show sig- 
nificant differences with respect to the data, refiecting large scale uncertainties. For specific choices 
of scales, the normalizations for either generator can be made to agree with the measurements. 
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The production of jets in association with vector 
bosons (y+ jets) in hadron colhsions is an important pro- 
cess in quantum chromodynamics (QCD) and is a sig- 
nificant source of background for many standard model 
measurements (e.g., ti production) and in searches for 
new phenomena (e.g., supersymmetry) . Such measure- 
ments at the Fermilab Tevatron colHder and the CERN 
Large Hadron Colhder (LHC) rely on accurate descrip- 
tions of F-|-jets production by particle- level event gen- 
erators. These models require validation with measure- 
ments of the properties of the F-l- jets system, especially 
as a function of jet multiplicity. 

In this Letter, we present new precision measure- 
ments of differential cross sections for the production of 
Z/j*+iets+X in pp collisions at a center-of-mass energy 
of 1.96 TeV. Cross sections are presented in bins of the 
transverse momentum {pt) of the N^^ jet in events con- 
taining at least = 1, 2, or 3 jets and are normalized 
to the measured inclusive Z/j* cross section to reduce 
uncertainties. The A'^ jets are ordered in terms of de- 
creasing Pt, and the Z/'y* is selected via its decay into 
an electron-positron pair. The data set corresponds to 
an integrated luminosity of 1.04 ± 0.06 fb~^ [l[. 

Previous Z/7*-|-jets measurements have focused pri- 
marily on measuring the jet multiplicities for up to 
three 0] and four [J jets. In addition, differential dis- 
tributions have been presented for the highest-pT (lead- 
ing) jet 0] and for the two leading jets [1]. In this Let- 
ter, we extend the measurements of differential distribu- 
tions by including the third leading jet, and by including 
a larger pi^* range. A major focus of the Letter is a 
comparison of the differential pT distributions to lead- 
ing order (LO) and next-to-lcading (NLO) pcrturbativc 
QCD (pQCD) predictions from MCFM as well as to 
results from several commonly-used event generators. In 
particular, we have investigated to what extent the dif- 
ferential Pt distributions for the higher jet multiplicities 
can be described by parton-shower-based event genera- 
tors like PYTHIA and herwig 0] , and how they com- 
pare to event generator predictions where matrix element 
and parton shower merging procedures are adopted, as 
in ALPGEN-I- PYTHIA and SHERPA @. For each pre- 
diction, the renormalization and factorization scale un- 
certainties are evaluated, and the choices of scales are 
modified to achieve an improved description of the mea- 
surements. The presented studies are of vital importance 
to understand the predictive power of the various event 
generator models for F-l-jets processes at both the Teva- 
tron and the LHC. 

The data set was recorded using the DO Run II de- 
tector, which is described in detail elsewhere Here 
we give a brief overview of the most relevant components 
for this analysis. The trajectory of charged particles are 
reconstructed using a silicon vertex tracker and a scin- 
tillating fiber tracker located inside a superconducting 
solenoidal coil that provides a magnetic field of approx- 



imately 2 T. The tracking volume is surrounded by a 
liquid-argon and uranium calorimeter, divided into elec- 
tromagnetic and hadronic sections with a granularity of 
An X A(f) = 0.1 X 0.1, where rj is the pseudorapidity 
[111 and (j) is the azimuthal angle. The third layer of 
the electromagnetic calorimeter has a finer granularity of 
A?7 X A(p = 0.05 X 0.05. The calorimeter consists of three 
sections, each housed in a separate cryostat, with a cen- 
tral section covering \r]\ < 1.1 and two end calorimeters 
extending the coverage to \ri\ « 4.2. Scintillators between 
the cryostats sample shower energy for 1.1 < \r]\ < 1.4. 

Electrons are identified based on their characteristic 
energy deposition signature in the calorimeter, including 
the transverse and longitudinal shower profiles. In addi- 
tion, a reconstructed track must point to the energy de- 
posit in the calorimeter, and the momentum of the track 
and the calorimeter energy must be consistent. Rejection 
against background from photons and jets is achieved 
with a likelihood discriminant which uses calorimeter and 
tracking information [l^ . 

Calorimeter jets are reconstructed using the DO Rmi II 
iterative seed-based cone jet algorithm [l3], using a 
splitting/merging fraction of 0.5 and a cone radius 
7^ ~ J(Ay)^ + (A(^)^ = 0.5, with y being the rapid- 
ity [1J|. The input objects are clusters of energy de- 
posited in the calorimeter. Rejection of jets arising from 
electronic noise in the calorimeter is achieved by us- 
ing quality and jet shape cuts. The reconstructed jets 
are corrected for the calorimeter response, instrumen- 
tal out-of-cone showering effects, and additional energy 
from multiple pp interactions and previous beam cross- 
ings. These corrections were derived by exploiting the 
Pt balance in 7-l-jet and dijet events. The measurements 
are performed for jets with pT > 20 GeV and |?7| < 2.5. 

Events were required to pass single or dielectron trig- 
ger requirements, and to contain two electron candidates 
with opposite sign electric charge, pt > 25 GeV, \r]\ < 1.1 
or 1.5 < I77I < 2.5, and a dielectron mass (Mee) satisfy- 
ing 65 < Mee < 115 GeV. A total of 65,759 events pass 
the selection before background subtraction. Of these, 
8,452/1,233/167 events have 1/2/3 jets or more, with 
pip" above 20 GeV. The efficiency for the trigger require- 
ments to be satisfied by Z/j*{—^ e+e~) events which ful- 
fill the other selection criteria was found to be ~ 100%, 
independent of the number of jets in the event. 

Backgrounds arising from events which contain two 
real electrons (e.g., WW, tt, and Z/^*{-^ r+r")) were 
estimated using event samples generated with pythia 
v6.323, with the underlying event model configured using 
Tune A [l5l. The events were passed through a GEANT- 
based [1^1 simulation of the detector response, and each 
event sample was normalized to higher-order theoretical 
predictions 0, [13 before being subtracted from the mea- 
sured data. For the inclusive data sample, the estimated 
sum of backgrounds arising from events containing two 
real electrons is ~ 0.2%. Of the background sources. 
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only ti production has an average jet multiplicity which 
is significantly larger than that of the signal process. As 
a result, ti is the dominant background for events con- 
taining two or more jets, contributing up to 6% (3%) to 
the measured data for large values of pT of the second 
(third) jet. The sum of all backgrounds containing two 
real electrons is estimated to be less than one third of the 
statistical uncertainty of the measured data in all bins of 
the measurement. By studying a sample of same-charge 
dielectron events, backgrounds arising from events with 
one or more mis-reconstructed electrons (e.g., W^-|-jets or 
multijet events) were found to be below 1% in the in- 
clusive Z/j*{-^ e~^e~) sample. For the signal samples 
containing at least 1,2, or 3 jets, no statistically signif- 
icant contribution from W-|-jets or multijet events was 
observed. For VF-|-jets, this was confirmed independently 
by using an event sample generated with pythia, yield- 
ing a contribution to the signal samples at the 0.1% level. 

The corrections of the reconstructed p!^^ spectra to the 
particle level [isj were determined using an event sam- 
ple generated with ALPGEN v2.05 + pythia v6.325 us- 
ing Tune A. The events were passed through a GEANT- 
based simulation of the detector response. The simulated 
events were overlaid with data events from random bunch 
crossings to reproduce the effects of detector noise and 
additional pp interactions. Particle-level jets were defined 
through the DO Run II iterative seed-based cone jet al- 
gorithm with TZ = 0.5 using a splitting/merging fraction 
of 0.5 and were required to satisfy \y\ < 2.5. The input 
objects were all stable particles except the two electrons 
defining the dielectron system and any photons in a cone 
oiTZ = 0.2 around the two electrons. The latter require- 
ment excludes most photons associated with the Z/^* 
decays, described in the simulation as QED final-state 
radiation. 

In the simulated sample, the ratio of the reconstructed 
p'^^ spectrum to the particle level spectrum defines the 
product of the efficiency of the detector (csim) and its 
geometrical acceptance (Agim)- The product (e x A) 
also includes the impact of the detector resolution on 
the reconstructed spectrum. To achieve agreement be- 
tween (e X A)sm-i and the corresponding quantity in data, 
(e X yl)data, the simulated event sample was modified in 
two steps. First, the simulated event sample was cor- 
rected so that its object identification efficiencies, energy 
scales, and energy resolutions correspond to those mea- 
sured in data. Next, the shapes of the p'^^ spectra in the 
simulated sample were reweighted at the particle level 
so that agreement with data was obtained at the recon- 
structed level. These two steps are explained in more 
detail below. 

The single electron identification efficiency was mea- 
sured in data for inclusive Z/^*{-^ e'^e~)+X events [l9l |. 
giving an efficiency, including statistical uncertainty, of 
^data = 0.77 ±0.01. The electron identification efficiency 
in the simulated sample, Cgj^j, was found to be w 10% 



higher than in data. To adjust the simulation to cor- 
respond to data, each reconstructed electron candidate 
in the simulated sample was rejected with a probability 
given by e^g^jj^/ e^jj^^. The electron efficiencies were binned 
in rje and As a cross check, the analysis was also per- 
formed using electron efficiencies binned in r^e and PT,e, 
resulting in changes of less than 1% in the final measure- 
ments. The same correction procedure was applied for 
the identification efficiency for jets, which was found to 
be ef^^^ = 0.98 ± 0.02 in data. 

The electron identification efficiency depends on the 
number of jets in an event and on their kinematics. Elec- 
tron candidates with a nearby jet are less likely to pass 
the electron identification criteria. To determine if the 
correlation between the electron efficiencies and the jet 
activity is correctly described by the simulation, two sets 
of comparisons were performed. First, the minimal AT?, 
separation between each electron and all jets in A^-jct 
events (A^ = 0, 1, 2, 3) in the simulated sample were com- 
pared with data. Second, the amount of soft QCD ra- 
diation not clustered into jets was studied using tracks 
which were not associated with any reconstructed elec- 
tron or jet. The multiplicity and pt sum of all such 
tracks in A^-jet events (A^ = 1, 2, 3), relative to those in 0- 
jet events, were compared between simulation and data. 
For both sets of comparisons, reasonable agreement was 
found. The total uncertainty due to differences between 
the simulation and the data in the dependence of the elec- 
tron identification efficiency on jet activity was estimated 
to be below 2% for all cross section measurements. 

Jet energy scale (JES) corrections were derived using 
7-|-jet and dijet events, assuring that the reconstructed 
p^^' in 7-l-jet events on average is equal to the particle- 
level p'^^ . The difference in calorimeter response between 
quark- and gluon-initiated jets introduces a dependence 
of the JES on the physics process. As a result, the JES 
corrections derived using 7-l-jet events do not guarantee a 
correct JES when applied to Z/7*4-jet events. However, 
the correction of data to the particle level depends only 
on data and simulation having a common JES. The factor 
(e X A) accounts for any differences between this common 
JES and the correct JES since particle-level jets in the 
simulated sample by definition have the proper energy 
scale. The px balance in Z/7*-|-jet events was used to 
adjust the JES in simulation to be equal to the JES in 
data. The correction is ~ 5% below 40 GeV and becomes 
negligible above 80 GeV. 

The jet energy resolution (JER) in data and simulated 
events was determined using Z/^*-\-]et events, and the 
reconstructed jet energies in the simulated sample were 
corrected to account for any differences. The JER dis- 
torts the steeply falling jet pT spectra, resulting in a net 
migration towards higher values of p'^^ . This leads to 
the reconstructed p^^' spectra being (5-15)% higher than 
they would have been for a detector with perfect jet en- 
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ergy resolution. The factor {ex A) accounts for this effect, 
and the measurements are fully corrected for the JER. 

After adjusting the simulated detector performance to 
the real detector performance, we tuned the p^* spectra 
in the simulation to those observed in data. For the mea- 
surement of each observable, the simulated event sam- 
ple was reweighted as a function of that observable at 
the particle level to ensure agreement of the distribution 
of the observable at the reconstructed level. After the 
reweighting, (e x A)sini is equal to (e x A)data within the 
uncertainties of the corrections applied to the simulated 
sample. The particle level spectra in data equal the re- 
constructed spectra in data divided by (e x A)sii„. The 



total efficiency for reconstructing a Z/j*{- 



e ' e 



jets-t-X event varies with pj, and is in the range 0.5-0.7 
for A^ 1 and 0.4-0.5 for N = 2,3. 

For events containing more than one jet, the A^*'' jet 
at the particle level might not be equal to the A^**^ jet at 
the reconstructed level. This can occur for events con- 
taining two jets with similar values of pi^'. The impact 
of this effect is part of (e x A) as defined above, and 
it is correctly taken into account if the JER is identi- 
cal for simulation and data, and if two jets have similar 
pip" values equally often in simulation and in data. The 
former is assured by the JER corrections applied to the 
simulated event sample, as described above. To test if 
the latter is satisfied, the ratios of the reconstructed pr- 
spectrum of the A^*'' jet to the spectra of the (A^ — 1)*'^ 
jet and the (A^ -I- 1)"^ jet in the simulated event sample 
was compared with those observed in data. Agreement 
was observed within statistical uncertainties. 

The final measurements are presented with two differ- 
ent particle-level selections in the mass range 65 < Mee < 
115 GeV: first without any further selections on the elec- 
trons (selection a) and secondly requiring p^ > 25 GeV 
and ly*^! < 1.1 or 1.5 < \y'^\ < 2.5 (selection 6). For each 
particle-level electron selection, the p!^ spectra were nor- 
malized to the inclusive Z/j*{—^ e^e^)+X cross section 
measured with the same particle-level selection. Selec- 
tion b corresponds to the kinematic range which is mea- 
sured in data. Selection a includes an extrapolation to 
the full range of lepton kinematics in order to simplify 
direct comparisons with other measurements. The ex- 
trapolation factor was derived from event samples gen- 
erated using SHERPA vl.1.1, ALPGEN v2.13 + PYTHIA 
v6.325 using Tune A and PYTHIA v6.323 using Tune A. 
The central value of the extrapolation factor was taken 
to be the SHERPA prediction, with the maximal devia- 
tion to the two other predictions being assigned as a sys- 
tematic uncertainty. The uncertainties due to the PDFs 
were evaluated using the Hessian method [l^l and were 
found to be negligible. The extrapolation from selection 
b to selection a increases the normalized differential cross 
section for the leading jet by w 10% below 100 GeV and 
decreases it by w 25% above 200 GeV. For the second 



(third) jet, the extrapolation changes the observable by 
less than 10% (2%). 

The systematic uncertainties of the measurements arise 
from the uncertainties of the background estimates, 
which were found to be negligible, and from the uncer- 
tainties of the corrections applied to the simulated event 
sample to assure that (e x A)su^ = (e x A)data- Each 
correction was varied separately within its uncertainties, 
and the resulting variations in the measured p^' spec- 
tra were added in quadrature to give the total system- 
atic uncertainty of the measurements. The largest source 
of uncertainty is the correction of the JES in the simu- 
lated event sample to correspond to the JES of the data 
sample, contributing (50-80)% of the total systematic 
uncertainty of the measurements. Additional uncertain- 
ties arise from the reweighting function applied to the 
simulated event sample at the particle level to ensure 
agreement with data at the reconstructed level, from the 
correction of the jet energy resolution and of the jet and 
electron identification efficiencies in the simulated event 
sample, and from the extrapolation from selection b to 
selection a. Presenting the measurements as ratios to the 
inclusive Z/7*(— > e^e~) + X cross section cancels the de- 
pendence on the uncertainty in the integrated luminosity 
of the data set. Additionally, most of the dependence on 
the uncertainties in the electron trigger and identification 
efficiencies also cancels. 

The choice of binning for the measurements was guided 
by the finite JER, which causes events to contribute to 
different bins at the particle and reconstruction levels. 
The purity of a bin is defined as the fraction of the sim- 
ulated events which are reconstructed in the same bin 
in which they were generated at the particle level. The 
widths of the measurement bins were chosen so that each 
bin has a purity of about 60%. 

The measurements presented above are compared with 
the predictions of several different theoretical models. 
For each model, the predicted jet px spectra are nor- 
mafized to the predicted inclusive Z/7*(— > e+e^) + X 
cross section. All predictions were generated using the 
CTEQ 6.1M [m parton density functions (PDFs) and 
the two-loop formula for the evolution of the strong cou- 
pling constant (as)- For the first and second jets the 
NLO MCFM predictions have been taken as the reference 
prediction; for the third jet, the leading-order (LO) MCFM 
prediction plays this role. The measurements and all the- 
oretical predictions are presented as ratios with respect 
to the reference prediction. 

The fully corrected measurements are summarized in 
Tables Hl lIIII and graphically represented in Figs. [TH31 
The data points are shown with statistical uncertainties 
(inner uncertainty bars) as well as with statistical and 
systematic uncertainties combined in quadrature (outer 
bars). Each data point is placed at the px value where 
the theoretical differential cross section is equal to the 
average cross section within the bin p^ . 
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TABLE I: The measurements of ^ ^ ^ x for the first jet (pr ordered) in Z/7*(— > e^e ) events with one or more jets. 



p^T > 25 GeV, ly"] < 1.1 or 1.5 < < 2.5, 

65 < Alee < 115 GeV 65 < M^e < 115 GeV 



Pt bin 


Bin center 


l/(7z/j* X dcr/dpy* 


5stat 


^sys 


l/az/^- X da/dp'^^ 


Sstat 


^sys 


[GeV] 


[GeV] 


[1/GeV] 


(%) 


(%) 


[1/GeV] 


(%) 


(%) 


20-28 


23.7 


6.81 X 10"'' 


1.6 


5.9 


7.19 X 10"^ 


1.6 


7.0 


28-40 


33.5 


2.99 X 10"^ 


2.1 


4.0 


3.22 X 10"^ 


2.1 


4.6 


40-54 


46.4 


1.23 X 10"^ 


3.1 


3.3 


1.37 X 10"^ 


3.1 


3.9 


54-73 


62.7 


5.04 X 10"" 


4.0 


3.3 


5.74 X 10"" 


4.0 


4.1 


73-95 


83.1 


2.03 X lO""* 


5.6 


7.6 


2.27 X 10"" 


5.6 


8.2 


95-120 


106.2 


7.29 X 10"^ 


8.4 


8.4 


7.62 X 10"^ 


8.4 


8.9 


120-154 


135.3 


2.64 X lO"'^ 


12 


10 


2.54 X lO"'^ 


12 


11 


154-200 


172.9 


8.08 X 10"'' 


17 


14 


6.99 X 10"'' 


17 


14 


200-300 


236.9 


7.46 X 10"^ 


42 


25 


5.58 X 10"^ 


42 


25 



TABLE IL The measurements of ^ ^ ^ x -j^ for the second jet [px ordered) in Z/'y*{—>- e+e ) events with two or more jets. 



Pt > 25 GeV, ly^^l < 1.1 or 1.5 < < 2.5, 

65 < Mee < 115 GeV 65 < M^e < 115 GeV 



Pt bin 


Bin center 


l/crz/j* X da/dp!^^ 


Sstat 


^sys 


l/fz/7* X da/dp^^^ 






[GeV] 


[GeV] 


[1/GeV] 


(%) 


(%) 


[1/GeV] 


(%) 


(%) 


20-28 


23.6 


1.30 X 10"^ 


3.7 


10 


1.39 X 10"^ 


3.7 


11 


28-40 


33.3 


4.23 X 10"" 


5.7 


5.2 


4.51 X 10"" 


5.7 


6.0 


40-54 


46.2 


1.57 X 10"" 


9.1 


6.2 


1.62 X 10"" 


9.1 


6.8 


54-73 


62.3 


4.17 X lO"'^ 


15 


8.5 


4.20 X lO"'^ 


15 


9.0 


73-200 


112.9 


2.96 X 10"® 


22 


7.4 


2.82 X 10"® 


22 


8.0 



TABLE III: The measurements of ^ ^ x -j^ for the third jet (pr ordered) in Z/^* {-^ e^e ) events with three or more jets. 







p^T > 25 GeV, \y^\ 


< 1.1 or 1.5 < 1? 


n < 2.5, 












65 < Mee < 115 GeV 




65 < Mee < 115 GeV 




Pt bin 


Bin center 


l/o"z/7* X da/dp'^^ 


5stat 


5sys 


l/f"z/7* X da/dp'^^ 


Sstat 


^sys 


[GeV] 


[GeV] 


[1/GeV] 


(%) 


(%) 


[1/GeV] 


(%) 


(%) 


20-28 


23.6 


2.22 X 10"" 


9.1 


14 


2.33 X 10"" 


9.1 


16 


28-44 


34.6 


4.40 X 10"^ 


17 


8.4 


4.48 X 10"^ 


17 


11 


44-60 


50.9 


8.67 X 10"® 


42 


11 


8.60 X 10"® 


42 


13 



The first comparison was performed between data and 
the pQCD predictions from mcfm v5.3 at NLO for the 
two leading jets, and at LO for the three leading jets. The 
central predictions were defined using factorization and 
renormalization scales fxp = I^^r ~ \J ^^z ^ Pt with 
Mz and pr.z denoting the mass and transverse momen- 
tum of the Z/7* boson. The sensitivity of the predicted 
cross sections to the choice of and /i/j was tested by 
varying their values up and down from the nominal value 
by a factor of two. The mcfm predictions were multiplied 
by correction factors accounting for multiple parton in- 
teractions (Cmpi) and hadronization (Cuad) before being 
compared to the measurements. 

The correction factors Cmpi and Cuad were estimated 
using inclusive Z/7*(— > e~'"e") event samples generated 



with PYTHIA v6.416 using Tune QW [23[, pythia v6.416 
using Tune SO HERWIG v6.510 + jimmy v4.31 [ist . 
ALPGEN v2.13 -I- PYTHIA v6.325 using Tune QW, and 
SHERPA vl.1.1. The central values quoted in Tables ITVl- 
IVII correspond to the predictions of pythia Tune QW. 
The maximal upwards and downwards differences be- 
tween PYTHIA Tune QW and the other four models are 
quoted as systematic uncertainties. 

Both the NLO and the LO mcfm predictions are in 
agreement with the measurements within the experimen- 
tal and theoretical uncertainties (Figs. [T]-[3|). At NLO, 
varying the scales up (down) by a factor of two changes 
the normalized 77)^* spectrum down (up) by factor of sal. 1 
for the leading jet, compared to a factor sal.2 for LO. 
For the second jet, the factors are wl.l (NLO) and wl.4 
(LO), and for the third sal. 6 (LO). These numbers illus- 
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FIG. 1; (a) The measured distribution of x ^^J^^^^j for the leading jet in ^/7*+jet+X events, compared to the predictions 

of MCFM NLO. The ratios of data and theory predictions to mcfm nlo are shown (b) for pQCD predictions corrected to the 
particle level, (c) for three parton-shower event generator models, and (d) for two event generators matching matrix-elements 
to a parton shower. The scale uncertainties were evaluated by varying the factorization and renormalization scales by a factor 
of two. 



TABLE IV: Correction factors for multiple parton interac- TABLE V: Correction factors for multiple parton interactions 



tions (Cmpi) and hadronization (Cnad) for 



'z/;- <ipr(l jet) • 



Pt bin 
[GeV] 



Cmpi ± (stat) ± (sys) Cnad ± (stat) ± (sys) 



20 - 28 
28 - 40 
40 - 54 
54 - 73 
73 - 95 
95 - 120 
120 - 154 
154 - 200 
200 - 300 



1.08 ±0 
1.04 ± 
1.02 ±0 
1.02 ±0 

1.01 ±0 

1.02 ±0 
1.04 ± 

1.03 ±0 
1.01 ±0 



UU_o.04 

LOOiHi 
m+0.00 

•"^-0.02 
ni+0.03 

nn+0.00 
•"^-0.03 
nq+0.00 

•"■^-o.o? 

nc;+0.02 

•"•J-o.oe 

•"^-0.05 



0.89 ±0 
0.90 ± 
0.90 ± 
0.92 ± 
0.93 ±0 
0.91 ±0 
0.92 ±0 
0.91 ±0 
0.92 ± 



'•""-0.03 

•""-0.01 
00+0.02 
•""-0.00 

•"^-0.03 
ni+0.01 
•"^-0.02 
nn+O.OS 
•"^-0.00 
no+OOS 
•"^-0.03 
O-+0.04 
•"*-0.06 
OO+0.05 
•"^-0.06 



trate the improved predictive power of the NLO compu- 
tation as compared with the LO one. The uncertainties 
of the MCFM predictions due to the PDFs were evaluated 
using the Hessian method. For the two leading jets, they 
vary from 5% at low to 10% at high pt^, and for the 
third jet they are found to be (5-15)%. 

Next, we compare the predictions of pythia v6.416 



(Cmpi) and hadronization (Cnad) for 



'Zl;-' dpT(2°'l jet) ■ 



Pt bin 
[GeV] 



Cmpi ± (stat) ± (sys) Cnad ± (stat) ± (sys) 



20 - 28 
28 - 40 
40 - 54 
54 - 73 
73 -200 



TTslTaoF^ 



+0.00 
-0.07 
-0.00 



1.10 ±0.01 
1.07 ±0.02: 

1.04 ± 0.03li^i!? 

1.05 ± 0.05t?, n2 



0.81 ±0.01 
0.83 ±0.01 



+U.OV 
-0.00 
+0.05 



-0.00 

0.85 ± O.Ollg-go 
0.87±0.03lg:oI 
0.83 ± 0.04+?,-,i,^ 



and HERWIG v6.510 + jimmy v4.31 with the measure- 
ments. These event generators describe jets through a 
parton shower using the approximation that parton emis- 
sions are soft or coUinear. For the hard qq —* Z/^* 
scattering, both generators use = Mz- For the 
parton shower, theoretical arguments favor the choice 
= a X with being the relative transverse 
momentum between the daughter partons in each 1^2 
parton splitting [2^ . This choice of fiR is adopted in both 
HERWIG and PYTHIA, with a = 1.0 being used in HERWIG 
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FIG. 2: (a) The measured distribution of 



for the second jet in Z/j* +2 jets+X events, compared to the 



predictions of MCFM NLO. The ratios of data and theory predictions to MCFM NLO are shown (b) for pQCD predictions 
corrected to the particle level, (c) for three parton-shower event generator models, and (d) for two event generators matching 
matrix-elements to a parton shower. The scale uncertainties were evaluated by varying the factorization and renormalization 
scales by a factor of two. 



TABLE VI: Correction factors for multiple parton interac- 



tions (Cmpi) and hadronization (Cuad) for 



■ dpT(3"^ jet) • 



Pt bin 
[GeV] 



Cmpi ± (stat) ± (sys) Cnad ± (stat) ± (sys) 



20 - 28 
28 - 44 
44 - 60 



1.15±0.02±H7 
1.10 ±0.031^4 

i.iiio.ioiHo 



0.76 ± 0.01+^;o^ 
0.81 ± O.OSlHo 
0.74 ± Omtom 



and in the final-state shower of pythia. For the initial- 
state shower, pythia using Tunc QW (Tune SO) sets 
a = (l-O)- Both HERWIG and pythia rewcight the 

leading parton shower emission to reproduce Z/7*+jct 
LO matrix-element computations [2^ . For the leading 
jet, PYTHIA using Tune QW shows a more steeply falling 
spectrum than observed in data (Fig. [1]). Tlie predic- 
tion of HERWIG + JIMMY shows good agreement with 
data at lowp^^*, but resembles pythia Tune QW at high 
p!^^. The change of slope around p^^* = 50 GeV can be 
traced back to the matrix-element correction algorithm 
in HERWIG [2^ . Comparisons to the measurements of the 



sub- leading jets (Figs.[2H3]) show that pythia using Tune 
QW and herwig predict more steeply falling p^' spec- 
tra than observed in data, in agreement with expecta- 
tions based on the limited validity of the soft/coUinear 
approximation of the parton shower. A newer pythia 
model with a pr-ordered parton shower, using Tune SO, 
gives a good description of the leading jet, but shows no 
improvement for the second or third jet. For the two 
pythia models, samples were generated with jip and /.i^ 
being varied up and down from the nominal value by a 
factor of two. As expected, decreasing and 11 r in- 
creases the predicted amount of events with one or more 
jets. The slopes of the predicted distributions do not 
change significantly as the scales are varied. 

Finally, we show comparisons with the ALPGEN v2.13 
-I- PYTHIA v6.325 and SHERPA vl.1.1 event generators. 
Both generators combine tree-level matrix elements with 
parton showers pol . [sol . 31 1, thereby utilizing matrix 
elements also for sub-leading jets. For the central 
ALPGEN+PYTHIA prediction, the factorization scale is 

given hy = \J M'^ + ^, whereas the renormaliza- 
tion scale is defined individually for each parton splitting 
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FIG. 3: (a) The measured distribution of x ^^^^^^..^ for the third jet in Zj'y* +3 jets+A' events, compared to the predictions 

of MCFM LO. The ratios of data and theory predictions to mcfm nlo are shown (b) for pQCD predictions corrected to the 
particle level, (c) for three parton-shower event generator models, and (d) for two event generators matching matrix-elements 
to a parton shower. The scale uncertainties were evaluated by varying the factorization and renormalization scales by a factor 
of two. 



using the CKKW prescription For SHERPA, both 

the factorization and the renormahzation scales are given 
by the CKKW prescription. For all three p^^* spectra, 
ALPGEN+PYTHIA predicts lower production rates than 
observed in data, but the shapes of the spectra are well 
described. SHERPA predicts a less steeply falling lead- 
ing 'p'^ spectrum than seen in data, leading to disagree- 
ments above 40 GeV. For the sub-leading p^* spectra, 
SHERPA predicts higher production rates than observed 
in data, but the shapes are well described. In agree- 
ment with Ref. [iOl, both ALPGEN-I-PYTHIA and SHERPA 
are found to show a sensitivity to the choice of scales 
which is similar to that of a LO pQCD prediction, re- 
flecting a limited predictive power. For the leading jet 
at pt ~ 100 GeV, the prediction of SHERPA with both 
scales shifted down by a factor of two is about three times 
higher than the ALPGEN-f PYTHIA prediction with both 
scales shifted up. This reflects both the size of the scale 
uncertainties and the difference in the central prediction 
between the two event generators. For ALPGEN+pythia, 
good and simultaneous agreement with data for all three 
leading jets is achieved through scaling /ij? and /ifl, down 



by a factor of two from the default values. For SHERPA, 
an improved description of data is achieved by scaling [ip 
and /i/j up by factor of two, but remaining disagreements 
with the measurements are seen for the leading jet below 
- 40 GeV. 

In summary, we have presented new measurements 
of differential cross sections for Z/"f*{-^ e+e~)-|-jets-|-X 
production in pp collisions at a center-of-mass energy of 
1.96 TeV using a data sample recorded by the DO de- 
tector corresponding to 1.04 ± 0.06 fb~^. The measure- 
ments are binned in the pt of the A^*'^ jet, using events 
containing at least N = 1, 2, or 3 jets, and are nor- 
malized to the measured inclusive Z/^*{-^ e+e") -I- X 
cross section. Predictions of mcfm at NLO, corrected 
to the particle level, are found to be in good agree- 
ment with data and have a significantly smaller scale 
uncertainty than MCFM at LO. The parton-shower based 
HERWIG and PYTHIA Tuue QW event generator models 
show significant disagreements with data which increase 
with ttJ^* and the number of jets in events. The newer 
PT-ordered shower model in pythia gives a good descrip- 
tion of the leading jet, but shows no improvement over 
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the old model for the sub-leading jets. The SHERPA and 
ALPGEN+PYTHIA generators show an improved descrip- 
tion of data as compared with the parton-shower-based 
generators. ALPGEN-|-pythia gives a good description 
of the shapes of the p!^^ spectra, while predicting lower 
production rates than observed in data. SHERPA pre- 
dicts higher production rates and a less steeply falling 
p'y^ spectrum for the leading jet than observed in data. 
For ALPGEN-j-PYTHiA, the factorization and renormaliza- 
tion scales can be chosen so that a good, simultaneous 
agreement with data is achieved for all three leading jets. 
For SHERPA, a similar level of agreement is achieved for 
the sub-leading jets, but some disagreements remain for 
the shape of the leading p^^* spectrum. Since the pre- 
sented measurements are fully corrected for instrumental 
effects, they can be used for testing and tuning of present 
and future event generator models. 
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